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Abstract 

The role of QCD scales and chiral symmetry in finite nuclei is examined. The 
Dirac-Hartree mean-field coupling constants of Nikolaus, Hoch, and Madland 
(NHM) are scaled in accordance with the QCD-based prescription of Manohar 
and Georgi. Whereas the nine empirically-based coupling constants of NHM 
span thirteen orders of magnitude, the scaled coupling constants are almost 
all natural, being dimensionless numbers of order one. We speculate that this 
result provides good evidence that QCD and chiral symmetry apply to finite 
nuclei. 
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Although QCD is widely believed to be the underlying theory of the strong interaction, 
a direct description of nuclear properties in terms of the natural degrees of freedom of that 
theory, quarks and gluons, has proven elusive. The problem is that at sufficiently low energy, 
the physical degrees of freedom of nuclei are nucleons and (intranuclear) pions. Nevertheless, 
QCD can be mapped onto the latter Hilbert space and the resulting effective field theory 
is capable in principle of providing a dynamical framework for nuclear calculations. This 
framework is usually called chiral perturbation theory (xPT). 

Two organizing principles govern this xPT: (1) (broken) chiral symmetry (which is man- 
ifest in QCD) and (2) an expansion in powers of (Q/A), where Q is a general intranuclear 
momentum or pion mass, and A is the generic QCD large-mass scale ~1 GeV, which in a 
loose sense indicates the transition region between the two alternative sets of degrees of free- 
dom indicated above (that is, quark-gluon versus nucleon-pion) . Typically, one constructs 
Lagrangians (that is, interactions) that display (broken) chiral symmetry and retains only 
those terms with exponents less than or equal to some fixed power of (1/A). The chiral 
symmetry itself provides a crucial constraint: a general term has the structure ~ (Q/A)^ 
and > is mandated. This guarantees that higher-order constructions in perturbation 
theory {viz., loops) will have even higher (not lower) powers of {Q/A). The price one pays for 
this mapping from natural to effective degrees of freedom is an infinite series of interaction 
terms, where coefficients are unknown and must be determined from experiment. 

To date only a few nuclear calculations have been performed within this framework. The 
seminal work of Weinberg highlighted the role of power counting and chiral symmetry 
in weakening N-body forces. That is, two-nucleon forces are stronger than three-nucleon 
forces, which are stronger than four-nucleon forces, • ■ • . This chain makes nuclear physics 
tractable. Van Kolck and collaborators developed a nuclear potential model, including 
one-loop (two-pion exchange) contributions. Friar and Coon developed non-adiabatic 
two-pion-exchange forces, while van Kolck, Friar and Goldman |Q examined isospin violation 
in the nuclear force. Rho, Park, and Min were the first to treat external electromagnetic 
and weak interactions with nuclei. Essentially all of this work was focused on few-nucleon 
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systems, where computational techniques are sophisticated. Only the work of Lynn on 
(nuclear) chiral liquids was specifically directed at heavier nuclei and, more recently, Gelmini 
and Ritzi ^ have calculated nuclear matter properties using lowest order nonlinear chiral 
eff'ective Lagrangians. 

Is there any evidence for chiral symmetry or QCD scales in finite nuclei? The tractability 
and astonishing success of the recent few-nucleon calculations of ^H, '^H, ^He, ^He, ^He, 
^He, ®Li, and ^Be with only a weak three-nucleon force and no four-nucleon force confirms 
Weinberg's power-counting prediction 0] and yields strong but indirect evidence for chiral 
symmetry. The work of Lynn established a procedure for going beyond few-nucleon 
systems. Nuclear (N-body) forces either have zero range or are generated by pion exchange. 
Following Manohar and Georgi [H we can scale a generic Lagrangian component as 
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where ip and vr are nucleon and pion fields, respectively, /^r and are the pion decay 
constant, 92.5 MeV, and pion mass, 139.6 MeV, respectively, A ~ 1 GeV has been discussed 
above, and {d^, m^r) signifies either a derivative or a power of the pion mass. Dirac matrices 
and isospin operators (we use t here rather than f) have been ignored. Chiral symmetry 
demands 

A = / + n- 2>0. (2) 

Thus the series contains only positive powers of (1/A). If the theory is natural @J^, the 
dimensionless coefficients cimn are of order (1). Thus, all information on scales ultimately 
resides in the cimn- If they are natural, scaling works. Our limited experience with nuclear- 
force models suggests that natural coefficients are the rule. 

Unfortunately, zero-range nuclear-force models are not widely used. However, a recent 
calculation has been performed using zero-range forces for an extended range of mass number 
A and this work provides significant new information on QCD and chiral symmetry in nuclei. 
Nikolaus, Hoch, and Madland (NHM) used a series of zero-range interactions to perform 
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Dirac-Hartree calculations in mean-field approximation for a total of fifty-seven nuclei. Their 
Lagrangian [using their notation] is given by 

— ^free ~l" -^4/ ~l" ^hot ~l" ^der ~l~ ^em i (3) 

where C/ree and £em are the kinetic and electromagnetic terms, respectively, and 



Chot = -l^sii^i^f - lisii^i^y - , and (5) 



Cder = -^55(5.^^) - ^5y(9,V'7^^)(5WV') • (6) 

In these equations, ip is the nucleon field, the subscripts S and V refer to the isoscalar-scalar 
and isoscalar-vector densities, respectively, and the subscripts TS and TV refer to the 
isovector-scalar and isovector-vector densities, respectively, containing the nucleon isospin 
operator r. The nine coupling constants of the NHM Lagrangian were determined in a 
self-consistent procedure that solved the model equations for several nuclei simultaneously 
in a nonlinear least-squares adjustment algorithm with respect to measured ground-state 



observables (Table IV of Ref. |10]). The predictive power of the extracted coupling constants 
is quite good both for other finite nuclei and for the properties of saturated nuclear matter 
(see Tables VIII, IX, and XI of Ref. |10|). 

£4/ contains four two-nucleon-force terms corresponding to A = 0, the first term of Chot 
is a three-nucleon-force term corresponding to A = 1, whereas the remaining two terms are 
four-nucleon-force terms corresponding to A = 2. Finally, C^er contains two nonlocal two- 
nucleon-force terms, also corresponding to A = 2. The derivative terms act on i/jiIj, rather 
than on just one of the fields, because the latter generate a factor E = M, the nucleon mass, 
whereas the former generate an energy difference that is considerably smaller. The latter 
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terms would spoil the series in Eq. (1) since M = A. However, either by a transformation 
or by rearranging the series, this problem could in principle be eliminated [^]. 

The construction of the NHM Lagrangian was motivated by empirically-based improve- 
ments to a Walecka type scalar- vector model , but using contact (zero-range) interac- 
tions to allow treatment of the Fock (exchange) terms. It was not motivated either by power 
counting or by chiral symmetry. The pion degrees of freedom are ignored and the Lagrangian 
is not complete; additional operators in each order of (1/A)^ are possible. Specifically, the 
NHM Lagrangian, Eqs. (4)-(6), has four operators in order (l/A)", one operator in order 
(1/A)^, and four operators in order (1/A)^, constituting an incomplete mix of three different 
orders in (1/A). 

Nevertheless, a meaningful comparison can be made of the generic chiral Lagrangian 
given by Eqs.(l) and (2) and the NHM Lagrangian given by Eqs.(4)-(6), precisely because 
our test of naturalness does not care whether a specific cimn coefficient is 0.5 or 2.0. Changing 
(refining) the model by adding terms would change all of the cimn, but the same test of 
naturalness still applies. Adding new terms would simply change a specific coefficient by an 
amount ~ 1 (or less). 

The nine coupling constants of the NHM Lagrangian are shown in Table 1, both in 
dimensional and dimensionless form [the latter obtained by equating Eqs.(l) and (4)-(6), 
with A = 1 GeV, using isospin operators t in Eq.(l), and solving for cimn in terms of a, /3, 
7, and 6]. In the former form they span more than thirteen orders of magnitude, while in 
the latter form six of the nine coupling constants can be regarded as natural. Only the very 
small aTs and large 75 and 'yv are unnatural. However, the sum of the latter appears to 
be natural, and we speculate that the difference may not be well determined in the least- 
squares adjustments to the measured observables. The unnaturally small axs, if correct, 
would presuppose a symmetry to preserve its small value. 

Although these results were not obtained as a test of chiral symmetry and QCD scales 
(NHM at that time were unaware of these developments) and hence are imperfect, they are 
conversely completely unbiased. This result is very indicative of the role of chiral symmetry 
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and QCD in finite nuclei. A systematic study of this approach is clearly indicated. 
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TABLES 

TABLE I. Optimized Coupling Constants for the NHM Lagrangian and Corresponding Di- 
mensional Power Counting Coefficients and Chiral Expansion Order 



Coup. Const. 


Magnitude 


Dimension 




Order 


as 


-4.508x10-^ 


MeV-2 


-1.93 


A° 


OiTS 


7.403x10"'^ 


MeV-2 


0.013 


A° 


ay 


3.427x10-^ 


MeV-2 


1.47 


A° 


OLTV 


3.257x10-5 


MeV-2 


0.56 


A° 




1. 110x10"^^ 


MeV-5 


0.27 


A-i 


IS 


5.735x10-1^ 


MeV-^ 


8.98 


A-2 


IV 


-4.389x10-1^ 


MeV-^ 


-6.87 


A-2 


Ss 


-4.239x10-1° 


MeV-^ 


-1.81 


A-2 


Sv 


-1.144x10-1° 


MeV-^ 


-0.49 


A-2 
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